Role of particle-number statistics in interference of independent Bose fields 
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We elucidate generally the interference of independent Bose fields in view of the conditional 
probability for the particle number measurements, and clarify its relation to the source number 
statistics. Despite lack of intrinsic phases, the interference phase can be inferred from the particle 
number registered at one detector by using the classical mean fields. If the conditional number 
distributions for the other detectors, given the outcome of the first detector, exhibit sufficiently 
narrow peaks around the values specified by the estimated phases, the mean field description is 
valid in a single run of interference. The widths in the conditional distribution are determined by 
the number statistics of the sources, among which notable scaling behavior is found depending on 
the detector configurations with the boundary at the Poissonian. The mean field description is found 
to be applicable to Poissonian and sub-Poissonian sources, whereas for super-Poissonian sources it 
is likely invalidated with the rather broad conditional distribution. 
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Interference is often considered as a signature of su- 
perposition in quantum systems. In particular, interfer- 
ence in many-body systems as a macroscopic quantum 
effect has been attracting many interests. In usual ex- 
periments, two states originating from a common source 
are subject to interfere, namely, each particle interferes 
with itself. However, in many-boson systems includ- 
ing lasers [1, 2] and atomic Bosc-Einstcin condensates 
(BECs) [3] , interference between independently prepared 
particles has also been observed. Such interference is of- 
ten explained with the spontaneous symmetry breaking 
for the relative phase, which gives nonvanishing expecta- 
tion values of the field operators or mean fields. In BECs, 
a U(l) symmetry is relevant for the global phase rotation 
of atomic wavefunctions, the breakdown of which relies 
on a nonphysical interaction [4, 5]. In optical systems, 
a U(l) symmetry also arises from lack of an absolute 
phase reference, which is ensured by the effective photon- 
number conservation in optical processes [6-8]. The U(l) 
symmetry breaking hence seems problematic in the ab- 
sence of real mechanism. 

The interference observed for independent sources un- 
der the U(l) symmetry has been attributed to the back- 
action of particle detection on the systems, which causes 
localization of the relative phase in a single run [6, 7, 9]. 
Another approach to the interference is to calculate the 
correlation functions of the particle numbers measured 
by the different detectors, which show the spatial mod- 
ulation. By evaluating the statistical moments of the 
Fourier components of the spatial modulation up to the 
fourth order, the plane- wave interference of atomic BECs 
is predicted in a single run with a random phase [4, 10]. 
This analysis exploits the nature of the plane- wave mode 
functions. Generally, some common understanding will 
be presented for the interference appearing under vari- 
ous configurations, which is based on the probability the- 
ory on quantum measurement. Moreover, there will be 
some intimate relationship between the interference and 



the particle-number statistics of sources, by considering 
the fact that the interference is observed so far for lasers 
(Poissonian states) and BECs (sub-Poissonian states). 

In this paper, we investigate the interference of in- 
dependent Bose fields under general configurations for 
sources and detectors, and clarify its relation to the 
source statistics. We examine the joint probabilities of 
the particle numbers registered by the detectors directly, 
rather than the correlations, to see the interference in a 
single run. The outcome at one detector provides infor- 
mation about the relative phase, despite lack of intrin- 
sic phases due to the U(l) symmetry. This information 
appears in the conditional distributions for the particle 
numbers at the other detectors, which are derived from 
the joint probabilities with the given outcome of the first 
detector. The relative phase is estimated by applying the 
mean field description to the measurement outcome. If 
the conditional distribution has sufiiciently narrow peaks 
around the values predicted by the estimated phases, it 
is almost certain that the outcome at the second detector 
takes a value close to one of the mean field predictions. 
Hence, the conditional distribution provides a quantita- 
tive criterion for the validity of the mean field descrip- 
tion. The mean field description is found to be appli- 
cable to Poissonian and sub-Poissonian sources, whereas 
for super-Poissonian sources it is likely invalidated with 
the rather broad conditional distribution. 

We consider a system of noninteracting Bose parti- 
cles, photons or cold atoms, where two independent 
sources are contained. The positive-frequency field oper- 
ator '(/'(x, t) is given generally in terms of the annihilation 
operators ai for a complete set of mode functions {^z}: 
^{ii.^t) = ai(j)i{ii.^t), where the time evolution of the 
free field is represented in the mode functions (l)i{x,t), 
which is determined in practice by expanding -0 alterna- 
tively in terms of the plane-wave modes. In order to de- 
scribe an interference experiment, the mode functions are 
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chosen suitably to provide the two independent sources as 
fli = d and 0,2 = b. For example, in interference between 
two wavepackets of light the wavevector distributions are 
localized around the central wavevectors of the respective 
sources. In the case of two atomic BECs [3], the initial 
mode functions (x, 0) are divided into two groups con- 
sisting of the eigenstates of the respective one-particle 
Hamiltonians with harmonic traps. In the following we 
assume for simplicity that all the particles are populated 
in the two source modes {I = 1,2), while the other modes 
{I > 3) are in the vacuum states. (This will be almost 
valid in typical interference experiments.) Then, the den- 
sity matrix for the sources is given hy p = pa® Pb, where 
each source state, respecting the U(l) symmetry, is given 
with the particle- number statistics Ps{N) [7] as 
00 

= 5]p,(7V)|iV)(iV| {s^a,h). (1) 

In the photon measurement for optical interference ex- 
periments, a commonly used photodetector records the 
number of photoelectrons emitted from the detector sur- 
face during a time interval T. The time and surface inte- 
grated photon-flux operator for the photoelectron emis- 
sion at the detector m is given [11] by 

Im=V'm dt dxdyij''{x,t)i'{x,t), (2) 

where rjm is the quantum efficiency, and the z axis is 
taken normal to the detector surface Sm- The band- 
width Ao; of the incident radiation is assumed to be 
small enough compared with the central frequency luq. 
The photon-flux operators in Eq. (2) are specifically ex- 
pressed as bilinear forms of the mode operators, /„i = 
i?|™''d|d;' , with the Hermitian matrices i?^™^ ob- 
tained from Eq. (2) by substitution -0^?/; 4)i4)ii . For the 
detection of cold atoms, we may take a resonant inter- 
action between the atomic internal levels and the probe 
light, which transfers the information of the atomic den- 
sity to ■ilj'^'tjj of the probe light [12]. Hence, the detection 
of cold atoms is treated in the same way as the photon 
number detection. 

The joint probabilities of the photon counts m, . . . , um 
by the M detectors (1 < M < Mend), which characterize 
the full statistics of interference, are given by 

P(ni,...,nM) = (:n ^am)"'"^-^"'" : V (3) 

where :: stands for normal ordering [11]. The flux oper- 
ators are presented explicitly as 

L = + Rt:^h% + R^h^b + 4:^6ta. (4) 

Here, it should be noted that the terms involving the 
vacuum modes {I > 3) are dropped in since they pro- 
vide null contributions to Eq. (3) as the normal-ordered 



expectation values. The mean particle number measured 
at each detector is given by 

{nrn) = (L) = Ri::^Na + Ri^'' Nb- (5) 

Here, Ng = Tr[jOsS^s] are the mean particle numbers ini- 
tially contained in the sources, which are assumed to be 
large enough to produce :§> 1 for high accuracy 

statistics. The coefficients Raa ^ and indicate the 

probabilities for each particle from the respective sources 
to fall into the detector m. They may represent the reso- 
lution of interference. Specifically, cx 1/Afond ^ 0, 
but keeping i?i™^./Vs 1 for (7i,„) :» 1, when the par- 
ticles are measured by almost continuously distributed 
many detectors, resulting in a fine interference pattern, 
e.g., spatial interference fringes [1, 3]. 

In the above sense, as seen in Eq. (5), a change of iii™' 
(or resolution) for the detectors may be viewed alterna- 
tively as an modification of the source statistics. Here, 
consider scaling of the detector matrices (by removing 
several detectors and changing the quantum efficiencies) , 

(g; M) = /q (m = 1, . . . , M) (6) 

with (g; M) = (m > M), and define the binomial 
distribution 

(0<A^'<iV). (7) 

In evaluating the joint probabilities, (:(/i)'^i • • • {ImY''' '■) 
contained in Eq. (3) are calculated for a Fock state 
\Na,Nb) with the normal-ordered expectation values 
((St)fc.(at)fcaa'=4'=') = [Na\/ {Na-kay]x[Nb\/ [Nb-Uy] 

{ka + kb = ki -\ h ^A/), which are multiplied by q'^-^q'^'' 

under the scaling. Then, by considering the relation 
q^iNl/iN - fc)!] = E^^'=fe Bf^,iq)[N'\/iN' - k)\], the ef- 
fccts of this scaling can be renormalized to the source 
statistics without changing the calculations in Eq. (3) as 

oo 

Ps{N;q)= E Ps{N')BS\q), (8) 
N'=N 

which is also normalized as the original pg{N). Hence, 
the number statistics of the sources may be replaced with 
the effective ones in Eq. (8) for any scaling of q, repro- 
ducing the same joint probability for the measurement 
by the AI detectors (namely the M-detector model): 

M),p,iN; q)} ^ P{n,, . . . , um). (9) 

This may be viewed as a renormalization transformation 
among the number statistics. It indicates universal re- 
lation for various interference phenomena, ranging from 
two-mode homodyne detection {M = 2) to measurement 
of spacial fringes (M = Mend ^1). According to Eq. (8), 
the mean Ns and variance Vg for the effective statistics 
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arc given in terms of the original ones as = qNg and 
Vs = q^Vs + {l — q)qNs. Then, for a sub- Poissonian distri- 
bution {Vs < Ng), the effective one is stih sub-Poissonian 

(V; < TV,) as 

VjNs^qiVjN,) + l-q. (10) 

The Poissonian form is preserved under the renormaliza- 
tion up to the scaling of mean as Ng — qNg. On the 
other hand, for a super-Poissonian distribution the effec- 
tive one is still super-Poissonian. 

We now examine the validity of the mean field descrip- 
tion for interference phenomena, where the field opera- 
tors are replaced with c-numbers as a — > a and b ^ (3 
(expectation values for coherent states |a,/3)). Specifi- 
cally, we have 

n,n = {a, P\Im\a, P) 

- (n™) + 2\R^^;:\nI/^nI'^ cosidab + Om), (11) 

where Na = ^6 = |/3|^, Sab = arga - arg^, 

9,n ~ argi?^™\ and (n,„) is the same as Eq. (5) for 
the U(l)-invariant sources. The set of {fim} exhibits the 
interference pattern with the cosine term in Eq. (11), 
which oscillates with 9„i depending on the detector loca- 
tion. The mean field description is, however, not directly 
applicable to the U(l)-invariant sources in Eq. (1) with 
(d^6) ~ 0, eliminating the cosine term in Eq. (11). Never- 
theless, by experiments and theoretical calculations the 
interference fringes are observed in a single run with a 
random relative phase for Poissonian sources (laser fields 
[1]) and sub-Poissonian sources (optical number states 
[6, 7] and BECs [3, 4, 10]). 

We hence consider the relationship between the in- 
terference phenomena and the source number statistics. 
Specifically, we examine the validity of the mean field 
description by inspecting the joint probability P{ni,n2) 
for any pair of detectors, say 1 and 2, depending on the 
source statistics. Given the outcome ni at detector 1, the 
mean field description in Eq. (11) provides an estimate 
for the relative phase, generally with two possibilities 5^^^ 
due to the cosine. Then, the outcome n2 at detector 2 is 
inferred with the estimated phases: 

ni=ni^ <5^b(rti) ^ n2[5^^{ni)]. (12) 

If the actual count n2 is close to one of n2[<5^b("'i)]' fixing 
the estimation of Job, we find that the interference occurs 
as described by the mean (classical) fields. This criterion 
for the interference can be checked readily by calculat- 
ing the conditional distribution Pc(n2\'ni) from P{ni, 712) 
with given ni. If Pc{n2\ni) has sufficiently narrow peaks 
at n2[d'^i,{ni)], the second outcome 712 should be close to 
either of the peaks with high probability. Specifically, 
the width of the peak should be no greater than that 
of the Poisson distribution e~"^(n2)"^/fi2!, which is the 



shot noise level for the coherent states |a,/3). Here, we 
conjecture that sub-Poissonian sources lead to the narrow 
peaks, showing the interference pattern. It is pointed out 
[13] that wavepackets emitted from a cavity maintain a 
pronounced relative phase coherence when the intracav- 
ity field has a narrow number distribution. Light beams 
from such sub-Poissonian cavities will exhibit the inter- 
ference. This phase coherence of each source is essen- 
tial to fix the interference phase in the number measure- 
ments. 

Consider first the case of fine detector resolution with 
i?'™'! <C 1 in the usual measurement of spatial inter- 
ference fringes. This case can be treated by scaling 
as the two-detector model with ^(^'-^^ = R^^'^'>/q ^ 1 
and g — > 0, which provides the same P[ni,n2) with 
the effective statistics in Eq. (9). Then, as seen in Eq. 
(10), the effective statistics of sub-Poissonian sources ap- 
proach the Poissonian for g ^ 0. Hence, by using any 
sub-Poissonian sources, essentially the same result is ob- 
tained for the interference fringes as the Poissonian case, 
where the mean field description is valid as numerically 
confirmed in the following. This is not the case for 
super-Poissonian sources. For i?^™^ = qw-™"^ — ^ with 
^(™) ^ 1 fixed, the large = R^/q 1/|-R(")|, which 
is required to produce {n.m) ^ 1, may derive even the 
larger V,, e.g., Vs oc Nf, for a super-Poissonian source, 
giving a nonzero q{Vs/Ns) for g ^> in Eq. (10). 

In order to examine the validity of the mean field 
description for general we have calculated nu- 

merically Pc{n2\ni) by using Eq. (3) for some typi- 
cal sources. The detector matrices are chosen for in- 
stance as = = 0.6i?, = i?ia' = OAR, 

^ib'^^l^ ~ Rad'^^ R'j,l'^\ giving the maximum interference 
term in Eq. (11), with the relative phase 6*2 — 6*1 = O.Ott. 
The first outcome is set as ni = 118, which corresponds 
to (m) +01 = Tl-39 and MSff^ini)] « 53,_113. 
Due to limitation on the numerical calculation, RNa = 
RNb — 100 arc taken, giving (ni) = (77,2) = 100 with 
Rad^^ + RH'^^ = R, and consistently (nj) + (712) = 200 « 
118 + (53 4- 113)/2. The scaling for the effective statis- 
tics is also used by taking R/q = R = 0.867 to calcu- 
late Pc("-2|'^i) for the increasing Na.b = 100/i? with the 
smaller R, after it is checked numerically for q ~ 0.5 
with Na^b = 25/i?. A bound on 82 — Oi may appear for 
the increasing R from the condition hi + h2 < Na + -^a 
{— 200/i?) due to the unitarity or the total number con- 
servation, e.g., 0.97r < 02 - 9i < n ioT R = 0.867. This 
is clearly seen in the familiar two-mode homodyne detec- 
tion, where /i,2 = (at ± &t)(a ± 6)/2 with e^C'^-^i) = -1. 

The resuhs for number states \N/R,N/R) with N = 
100 and some values of R are shown in Fig. 1. The 
case of Poissonian source is also plotted for comparison, 
corresponding to i? — 0, where the Poisson distribu- 
tion oc {h2)"^/n2l for 712 is confirmed around the peaks 
(though rather broad due to not so large (77.1.2) = 100). 
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FIG. 1. (Color online) Conditional distribution for number 
states |iV/-R, N/R) with iV = 100 and some values of R. The 
Poissonian case corresponds to _R — )■ 0. The mean field values 
?i2[5^^(ni)] ~ 53, 113 for rii — 118 are indicated with vertical 
dotted lines. 
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FIG. 2. (Color online) Conditional distribution for super- 
Poissonian sources with 'P{a) in Eq. (13). The mean field 
values are shown the same as in Fig. 1. 



The peaks agree with «2[i^^fc(^i)] ~ 53, 113 (vertical dot- 
ted lines), and exhibit the narrower widths than the Pois- 
sonian case. Therefore, the mean field description is valid 
for these sub-Poissonian number states and also their ef- 
fective statistics, i.e., the binomial distributions in Eq. 
(7). Here, the limit i? — > 1 becomes unphysical with the 
dominating n2[i5^ft(??i)] to give ni + n2 w 118 + 113 > 
200(i? = 1), violating the unitarity. 

We have also considered a super-Poissonian source 
with a U(l)-invariant 'P-representation as 



Via) oc {\af /QN^^^-^ eM-\af /QN) 



(13) 



where Q ^ {V - N) /N^ with Q > 0. The limit Q ^ 
corresponds to the Poissonian, whereas Q = 1 to the 
thermal state. The conditional distribution is shown in 



Fig. 2, which does not depend on R in this case with 
RNa.b (=100) fixed. The increasing Q broadens the dis- 
tribution, eventually washing out the peaks. We have 
further examined the single-photon-added thermal state 
[14]. This nonclassical super-Poissonian state has the 
variance smaller than the thermal case. Despite this fact, 
for the small R [RNa.b fixed), the conditional distribu- 
tion becomes flatter than that for the thermal sources. 
These results indicate that the behavior of interference 
is rather complicated for super-Poissonian sources, likely 
invalidating the mean field description. 

To conclude, in view of the conditional probability 
for the number measurements, we have elucidated the 
common mechanism for the interference of independent 
Bose fields under various situations, ranging from two- 
mode homodyne interference to spacial fringes. The in- 
terference is determined by the source number statistics, 
among which the scaling behavior is present depending 
on the detector characteristics with the boundary at the 
Poissonian. For sub-Poissonian and Poissonian sources 
the interference pattern appears in a single run, consis- 
tently with the mean field description, whereas this is not 
the case for super-Poissonian sources. It will be a chal- 
lenge for future experiments to confirm the role of source 
statistics with the scaling behavior, by preparing various 
source states and detector configurations. 
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